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Abstract. HD 553 was discovered as an eclipsing binary by Hipparcos. Here, we present the first radial velocity 
curve for this system. It is found, that HD 553 is a double-lined spectroscopic binary. Despite the large luminosity 
difference, the two components of this system are of very similar mass. The primary, a KO-giant, fills a large 
fraction of its Roche-lobe. The secondary is, despite its very similar mass, still a late-type dwarf. The radial 
velocity curve allows to constrain several stellar and system parameters. 
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1. Introduction 

Additional to the astrometric parameters of positions, 
proper motions and trigonometric parallaxes, the highly 
successful Hipparcos mission also provided photomet- 
ric information for the target stars. Thus, a large 
number of new variables has been discovered by 
Hipparcos. Among these is HD 553 (=HIP 834; 02000 = 
00 h 10 m 10!4647, 62000 = +64° 38' 48'.' 174; V« 8 m 2). 
Kazarovetset al. ( 1999| ) subsequently classified the system 
as an eclipsing binary of the (3 Lyr type (i.e. two minima 
are present and there is no constant part due to the ellip- 
ticity effect) and gave it the name V 741 Cas. Strassmeier 
et al. ( [2000D identified HD 553 as a single-lined spectro- 
scopic binary and included it into their list of potential 
targets for Doppler-imaging from their survey of stellar 
activity indicators. 

The Hipparcos catalogue (ESA [1997 ) provides the fol- 
lowing information about this system: 

- V=8 m ll, B-V=l m 032 ± m 016, V-I=l m 00 ± m 01 

- 7T = 4.74 ± 0.83 mas 

- period P = 9.0576 ± 0.0008 days, T o =2448504.0610 as 
time of zero phase (from fits to the Tycho-lightcurve). 

HD 553 is an ideal target for Doppler imaging. At ob- 
servatories north of latitude 26°, like the NAO Rozhen, 



Send offprint requests to: L. Iliev, liliev@libra.astro.bas.bg 

* Based on observations made with the Coude-spectrograph 
at the 2m-RCC-telescope of the National Astronomical 
Observatory Rozhen, Bulgaria. 



Bulgaria (ip rj 42°), it is circumpolar, allowing for in prin- 
ciple continuous observations. This would allow to follow 
the spot evolution on any timescale (longer than the rota- 
tion period) without the often 1 yr timegap between subse- 
quent maps for other Doppler imaging targets. Assuming 
that the orbital and rotational axes of the eclipsing binary 
are aligned, the inclination of the rotational axis is high 
allowing to observe most of the stellar surface. It is bright 
enough to allow for good signal-to-noise spectra within a 
reasonably short time. Because the eclipse, when the pri- 
mary is in the back, is far from total we can also resolve 
most of the stellar surface during the eclipse. Its projected 
rotational velocity v sin i = 42.5 kms -1 (see Sect. |^) is 
high enough to allow for reliable Doppler imaging even 
with a moderate resolution of about 30 000, as obtainable 
with the Coude spectrograph at Rozhen. Finally, the pe- 
riod of 9 days is long enough that one good spectrum per 
observing night is enough to provide a good phase cov- 
erage, and thus the object can easily be included in any 
spectroscopic night's object list. On the other hand, the 
period is short enough to allow the data for a complete 
map to be collected within the typical 1-2 week spectro- 
scopic observing runs. 

Therefore, we became interested in this object and 
started an observing campaign. The first step for reliable 
Doppler imaging is to find a good radial velocity curve. 
The problem of Doppler imaging is ill-posed and has a lot 
of free parameters. Thus, any additional constraints that 
can be put prior to Doppler imaging is of great value. Since 
Doppler imaging relies on distortions in the line profiles 
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caused by the surface structures, it is important to remove 
any other distortion of the line profiles, in particular any 
shift of the lines with respect to the same lines at other 
phases. In principle, the radial velocity shifts from the 
orbital motions could be removed prior to Doppler imag- 
ing by simply cross-correlating the spectra. However, the 
spots on the stellar surface introduce systematic errors in 
the shifts obtained in this way. It is therefore useful to 
determine the radial velocity based on a large time scale; 
each individual radial velocity will still be influenced by 
the line profile distortions, but the limited lifetime of the 
spots ensures that spectra at the same phases taken a 
long time apart average out these distortions. Also, any 
reliable determination of the orbital parameters is based 
on a good value for the orbital period (all other parame- 
ters more or less correlate with the orbital period) , which 
becomes more and more accurate the longer the time base 
of the radial velocity curve. 

Here, we report on the first radial velocity curve of 
the HD 553 system based on 20 high-accuracy radial ve- 
locities obtained over the time base of 2 years, and the 
conclusions about the system and stellar parameters that 
can be drawn from it. 



2. Observations and reductions 

The observing material consists of 20 spectra ob- 
tained with the Coude spectrograph at the 2m-RCC- 
telescope of the National Astronomical Observatory 
Rozhen, Bulgaria. The data were collected in 4 observ- 
ing runs in October/November 1999, April and October 
2000, and September/October 2001. 

The grating 632/22?3 and an entrance slit width of 300 
lira were used, resulting in a resolution of about 30 000. 
The spectra were recorded with the Photometries AT200 
camera on a 1024x1024 CCD with pixel size 24x24/^m 2 , 
yielding a spectral range of about 10 nm for each spec- 
trum. In different seasons, different spectral regions have 
been observed: in 1999, the region around the Cal line at 
643.9 nm, in 2000 and 2001, the region around the Fei line 
at 617.3 nm have been observed; additionally, in 2001, an- 
other region, that around a group of telluric lines at 630 
nm, was observed immediately before or after the Fe I- 
region, which allows an independent check of the wave- 
length scale via the atmospheric lines and of the radial 
velocities (RVs) via the comparison of the RVs from the 
different regions at virtually the same time. The average 
signal-to-noise ratio of the spectra is 136. 

Additional to the spectra of HD 553, several spectra 
of the RV standard (3 Gem have been obtained in each 
observing run, except for Oct 2000, when we observed the 
RV standard a Cas. Both stars are classified as K0III, and 
are therefore suitable templates for HD 553. 



The information for each spectrum is given in 
Table |J 

Each observed stellar spectrum is accompanied by a 
pair of comparison spectra (one taken before, the other 
after the stellar exposure) and a number of flatfield im- 
ages. In several nights, a large number of bias-images was 
obtained, which were found to be unchanging within an 
observing run; they were averaged, then subtracted from 
each image of the same run. 

The reduction of the material has been done using 
the spectrum reduction package 4A (Ilyin |2000 ), and 
consists of the steps bias image subtraction, flatfielding, 
scattered light subtraction, optimal spectrum extraction, 
wavelength calibration. 

The wavelength calibration is based on the comparison 
spectrum images consisting of Th-Ar emission lines. Each 
Th-Ar-image consists of a pair of comparison spectra: one 
is above, the other below the area covered by the stellar 
spectrum, but overlapping with it. This allows to take into 
account a slope of the lines of constant wavelength across 
the stellar spectrum. The observation of a comparison im- 
age before and another one after the stellar exposure al- 
lows to take into account a possible systematic drift of the 
pixel coordinate system along the wavelength axis. Thus, 
the final dispersion curve of each stellar image is based 
on a total of 4 Th-Ar-comparison spectra, and is interpo- 
lated to the central position of the stellar spectrum and 
the time of mid-exposure of the stellar spectrum. The dis- 
persion curve is approximated by a polynomial of degree 

3 (except for Sept/Oct 2001, where a degree of 5 was nec- 
essary) and has typically an accuracy of 0.02-0. 04 kins" 1 . 

After the wavelength calibration, the spectra are trans- 
formed into heliocentric wavelengths and normalized to 
the apparent continuum which is a polynomial of degree 

4 approximating the high points of the stellar spectrum. 

It is worth mentioning that 4A estimates the errors 
of each pixel in the original image (after bias subtrac- 
tion) assuming Poissonian noise and propagates these er- 
rors through each step of the reduction procedure. In this 
way, reliable errors are available for each pixel in the fi- 
nal spectrum in both flux and wavelength resulting also 
in individual errors for each RV. 



3. Radial velocity measurements 

Radial velocities were determined for each spectrum by 
cross-correlation with the average of all spectra of the ra- 
dial velocity standard observed during the same observ- 
ing run. For this, the spectra of the standards were ar- 
tificially spun up to different values of the projected ro- 
tational velocity v sin i, and the spectrum which max- 
imizes the cross-correlation peak height gives the best 
value of v sin i. This was done for both standards (3 Gem 
and a Cas and for all spectra of HD 553. The result is 



1 Table jj] is also available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/Abstract.html 
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Fig. 1. The spectrum of HD 553 (thick line, top) com- 
pared to the spun-up spectrum of /3 Gem (v sin i = 
42.5 km s -1 ; thin line), both shifted up by 0.5. At the 
bottom is the difference between these two spectra, 
showing some weak, narrow absorption lines of the sec- 
ondary. The central spectrum is the original spectrum of 
(3 Gem, used to measure the RV of the secondary, which 
is here -101 km s^ 1 with respect to the primary (HJD 
2452185.5178). The spectra of HD 553 are in the restframe 
of the primary (determined from a preliminary orbit fit, 
based only on the primary's RVs). 



v sin i = 42.5 ± 0.5 km s" 1 , independent of which stan- 
dard is used, independent of the spectral region and of 
the phase in the binary orbit. The latter allows the con- 
clusion that the ellipticity of the primary in the binary 
(the giant) is not very pronounced. 

Before cross-correlating the artificially broadened 
spectra of the templates with the spectra of HD 553 they 
were corrected for the radial velocity of the corresponding 
standard star. For (3 Gem, RV=+(3.3±0.1)kms _1 accord- 



ing to the Astr onom ical Almanach for 2002 ( 2001 ) and 
Stcfanik et al. ( 1999 ). For a Cas, in recent publications 
several different RVs are given, casting some doubt on 
its suitability as a radial velocity standard: RV=— (3.9 ± 
0.1) k ms" 1 according to the Astronomical Almanach for 
2002 ( |200 1[) RV=-(4.3 ± 0.2) kms" 1 according to Udry 
et al. ( [1999D , and RV=-( 4.79 ± 0.04) kms" 1 according 
to de Medeiros & Mayor (199E). Our investigation (see 



Sect. |]) shows that the value by de Medeiros & Mayor 
(1999) is best. Note, that by observing the standard stars 
in the same run as the spectra of HD 553 (usually im- 
mediately before or after HD 553) and reducing them in 
exactly the same way, then moving them to RV=0 using 
the values from the literature, the RVs of HD 553 are on 
the system of the literature velocities of the standards, 
without any contribution from some systematic velocity 
zero-points (which could be different in different runs) (see 
also Sect. |). 

We also searched for signs of the secondary in the spec- 
tra. Indeed, a faint spectrum displaying narrow lines mov- 



ing in antiphase to the strong, broad lines of the giant 
were found, after the spectrum of (3 Gem had been sub- 
tracted from that of HD 553. The RVs from these lines 
were measured by cross-correlation with the original (not 
spun-up) spectrum of the RV standard of that run. RVs 
of the secondary are not measureable when the secondary 
is eclipsed. Figure [l] shows the spectrum of HD 553 in the 
617.3 nm region and identifies some lines of the secondary. 
It can be seen that the continuum of the broadened (3 Gem 
spectrum tends to be below that of HD 553. The reason is 
the imperfect continuum normalization for HD 553, which 
is due to the very low signal-to-noise and to the fact that 
its (pseudo-) continuum is only defined by a very few re- 
gions without broad, blended features. However, the dif- 
ference spectrum is mostly flat so that RV measurements 
should not be influenced and no attempt has been made 
to correct for this. Note that some "emission" features 
appear in the difference spectrum which result from an 
imperfect match of some of the features in the artificially 
broadened (3 Gem spectrum and in the HD 553 spectrum. 
Yet, the lines of the secondary are clearly visible and mea- 
sureable. 

All the strong lines in the (3 Gem spectrum are repre- 
sented in the spectrum of the secondary; several weaker 
features are also identifiable. The lines are narrow, which 
can be seen e.g. with the clearly resolved double feature 
near 613.5 nm. This indicates a much lower v sin i of the 
secondary; together with the low luminosity and the pres- 
ence of the same strong lines as in the spectrum of the late 
giant (3 Gem, this makes the secondary a late-type dwarf. 
Therefore, we decided to use the unbroadened spectum 
of j3 Gem as a template for the RV measurements of the 
secondary. Earlier experience has shown that a good RV 
of (3 Gem can be measured with the Sun as a template. 
So, a mismatch of the spectral types of (3 Gem and of the 
secondary should have only a small effect on the final RV 
of the secondary. 

Unfortunately, the lack of proper standards, the com- 
paritively low signal-to-noise ratio in the spectrum of the 
secondary, and the also visible systematics resulting from 
the not perfect match of the broadened (3 Gem spectrum 
and the HD 553 spectrum make it impossible to analyse 
the spectrum of the secondary in more detail. 

All RVs are given in Table [l] together with their errors 
as determined by 4 A. These take into account the errors of 
the pixels in the spectra in the overlapped region of both 
the programme and the standard spectrum as well as the 
height and curvature of the cross-correlation maximum 
(for details, see Ilyin 2000). 



4. Orbital solution 

The radial velocities in Table ^ were used in a double-lined 
orbit fit. A few details about the fit routine can be found 
in Duemmler et al. (1997). The weights of the individual 



RVs were obtained as the inverse variances, based on the 
RV errors. However, it turned out, that the primary and 
secondary had very different values of % 2 ; therefore, the 
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Table 1. The radial velocities of HD 553. The columns 
are: The heliocentric Julian date (HJD-2450000) of the 
mid-exposure, the approximate centre of the spectral re- 
gion in nm, the signal-to-noise ratio (SNR) of the spec- 
trum, the radial velocities and errors of the primary and 
secondary, respectively, in kms -1 . 



HJD-2450000 


A c 


SNR 


RVi 


RV 2 


1478.6010 


644 


76 


17.26±0.38 


-97.45±0.77 


1481.3821 


644 


162 


-88.47±0.38 


11.87±0.66 


1482.4879 


644 


157 


-102.36±0.26 


25.36±0.69 


1483.3475 


644 


189 


-89.04±0.23 


12.81±0.55 


1484.3707 


644 


151 


-51.71±0.20 


— 


1485.3671 


644 


174 


-9.59±0.23 


-70.83±0.54 


1648.5779 


617 


73 


-3.36±0.40 


-74.37±0.69 


1652.3007 


617 


59 


-40.38±0.32 




1831.3489 


617 


160 


23.94±0.35 


-103.73±0.53 


1832.2459 


617 


121 


8.41±0.36 


-89.52±0.68 


1832.6138 


617 


182 


-3.91±0.34 


-77.18±0.50 


1833.3567 


617 


160 


-33.31±0.29 




1834.3547 


617 


169 


-74.61±0.36 


-0.42±0.53 


2181.4355 


630 


120 


-73.98±0.44 


-3.81±1.63 


2181.4664 


617 


125 


-71.79±0.26 


-3.30±0.54 


2182.5348 


630 


126 


-28.18±0.33 




2182.5675 


617 


147 


-26.87±0.19 




2183.5393 


630 


102 


10.27±0.45 


-89.74±1.18 


2185.4905 


630 


145 


11.20±0.36 


-90.93±1.06 


2185.5178 


617 


132 


10.69±0.27 


-90.48±0.37 



weights of the secondary were reduced by a common factor 
Xl/X% = 0-108. 

While the large ratio of the \ 2 values of the primary 
and secondary have shown that the errors for the sec- 
ondary given in Table ^ are clearly underestimated, the 
overall reduced \ 2 of the fit of 2.1 is quite good. Although 
the formal probability of exceeding a value of 2.0 is just 1% 
for 15 degrees of freedom, experience with many fits has 
shown that, even for larger numbers of RVs with smaller 
errors, values of 2-3 for \ 2 are not uncommon. However, 
larger deviation than expected from the errors may also re- 
sult from stellar activity: starspots distort the line profiles, 
and cause the RV to deviate systematically. Therefore, we 
have investigated the O-C values as a function of phase 
and observing run, but we do not find any systematic (see 
Fig. ||). One should, however, keep in mind that 20 mea- 
surements is still a rather small statistics. 

The errors of the fit-parameters were determined by 
10 000 bootstrap runs (see Duemmler et al. 1997 and, e.g., 
Efron & Tibshirani 1993). To each RV computed from the 



orbital solution we add a randomly chosen value taken 
from the O-C values. Since on average low weighted RVs 
tend to have larger O-C, each O-C is accompanied by its 
original weight; the choice of O-C values is done with re- 
placement, i.e. the same O-C can be chosen several times 
while others are not used. With this artificial data set a 
new orbital fit is performed, leading to a slightly different 
set of parameters than the fit to the original data. The 



Table 2. The orbital parameters of the double- lined spec- 
troscopic binary HD 553. Errors of the parameters are 
based on 10 000 bootstrap runs. The primary derived pa- 
rameter is P res t) the period in the restframe of the bi- 
nary system, i.e. the observed period corrected for the 
Doppler-effect caused by 7. All dynamical derived param- 
eters are based on P res t, while the epochs To, T con j are in 
the observer's system, i.e. based on P bs- To is the time 
of maximum RV for the primary (the giant), T con j is the 
conjugation time with the primary in the back. \ 2 is the 
reduced \ 2 °f the fit, i.e. the weighted sum of squared de- 
viations, divided by the number of degrees of freedom, a 
is the mean standard deviation of a single RV with aver- 
age weight, <j\ and 172 are the same for the primary and 
secondary, respectively. 



parameter 



K 2 (kms" 1 ) 
Pobs (days) 
7 (kms -1 ) 
T (HJD) 

Prest (days) 

q — 7712 /mi 
ai sin i (Rq) 
a 2 sin i (Rq) 
mi sin i (Mq) 
m 2 sin i (Mq) 
T conj (HJD) 



a (kms -1 ) 
a 1 (kms -1 ) 
o 2 ( kms -1 ) 



63.32±0.19 
66.04±0.41 
9.05997±0.00011 
-38.93±0.12 
2451813.0877±0.0036 

9.06115±0.00011 
0.9589±0.0066 
11.336±0.034 
11.822±0.073 
1.037±0.013 
0.995±0.009 
2451815.3527±0.0035 

2.14 
0.56 
0.59 
0.97 



error of a parameter is given by the standard deviation 
of this parameter from all the bootstrap samples. The re- 
sulting orbital and derived parameters are presented in 
Table |. The orbit is plotte d in Fig. |[ 

Strassmeier et al. ( 2000 ) have also measured two RVs 
for HD 553. While these two RVs increase the time span of 
the RV curve by another year, the large errors compared 
to our RVs lead to very small weights in the fit. The result 
is identical to the fit presented in Table 0, except that the 
errors of the parameters even increased slightly. We there- 
fore decided to ignore the Strassmeier et al. measurements. 

We also attempted to fit an eccentric orbit. The result 
e = 0.0016 ± 0.0023 with u = 350?7 ± 0?1 clearly shows 
that the or bit is circul ar. Th e Lucy-Sweeney P-test (Lucy 
& Sweeney 1971 , Lucy 1989 ) gives a 96.2 % probability for 
a circular orbit. 

Another fit allowed for different 7 velocities for the 
4 observing runs, while forcing all other parameters to be 
the same for the 4 data sets. That means that we are look- 
ing for systematic RV errors between the different runs, or 
for systematic changes of the 7 velocity from run to run, 
which could be caused by a third body. Within the errors, 
the orbital parameters do not change compared to those 
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HD 553: Radial velocity cu 




-D.2 0.0 



HD553; RV curve, Q-C 



I - 



Fig. 2. top: The radial velocity curves of the two stars 
in the HD 553 system. Filled symbols are used for the 
primary (the giant), open symbols for the secondary; cir- 
cles denote the data for Oct/Nov 99, diamonds those for 
Apr 00, triangles those for Oct 00, and squares those for 
Sep/Oct 01. Errorbars are typically smaller than the sym- 
bol size. 

bottom: O-C vs. phase. 

given in Table |2|; the results for 7 are: 



Oct/Nov 1999 71 = 

Apr 2000 72 = 

Oct 2000 73 = 

Sep/Oct 2001 74 = 



-39.27 ±0.19 
-38.62 ±0.34 
-38.73 ±0.25 
-38.75 ±0.26. 



(1) 



The largest difference, 71 — 72, is only 1.7 a, i.e. there 
is no indication for any significant change of 7 over the 
two-year interval or for a significant error in the veloc- 
ity zero points of the different runs. It is noteworthy 
that 73 is based on RVs measured with a Cas as stan- 
dard, for whose RV we adopted the value given by de 
Medeiros & Mayor ( 1999| ). Would we have adopted the 
other values given in the literature, 73 would deviate by 



0.5kms _1 (Udry et al. 1999), respectively; it is clear, that 
the value by de Medeiros & Mayor ( 1999 ) is best, at least 
at the time of our measurements, Oct 2000. 

It is remarkable that the orbital parameters in Table || 
show a system with almost identical masses for the two 
stars, although one of the stars is clearly evolved. 

It should be mentioned that our period P bs is longer 
than both the period given in the Hipparcos catalogue 
(ESA |1997] P = 9"?0576 ±0^0008) and the period 
given by Strassmeier et al. ( ^OOOfc P = 9 d 0580). While 



Strassmeier et al.'s period is consistent with the one given 
by Hipparcos, our period is longer by 3c. However, since 
the Hipparcos period is based on a curve fit, whose shape 
is not specified, we feel that our period is more reliable. 

5. Conclusions and discussion 

We have found the Hipparcos-detected eclipsing binary 
HD 553 to be a double-lined spectroscopic binary with a 
giant as primary and a faint companion of almost equal 
mass. From 20 high-quality RVs collected over the course 
of 2 years, the radial velocity curve was derived, leading to 
an improvement of the binary period and the derivation 
of other orbital and physical parameters. 

If we assume that the primary is in synchronised ro- 
tation (i.e. the system is old and the primary has been a 
giant for a long time), the rotation and orbital axes should 
also be aligned, i.e. P ro t = Porb, for which the period P res t 
in the rest of frame of the system should be used and 
*rot = *orb — Then, the projected radius is given by 
v sin i via 



Ri sin i 



Pv sin 
2tt 



(7.608 ±0.090) i? , 



(2) 



with v sin i = (42.5±0.5) kms -1 as determined in Sect. |[ 
The assumption of aligned axes and synchronisation of 
rotation and orbit could, however, be wrong. While it is 
true that the system is old for the primary to become a 
giant by now, so possibly the axes are aligned, the periods 
may be different. During the expansion of the primary to- 
wards the giant branch the rotation period should increase 
due to angular momentum conservation. A longer rotation 
period would increase the radius given above. It would 
also imply that the primary's expansion happened quite 
recently, so that the tidal forces were not able to speed up 
again the primary's rotation. This is highly likely, because 
the secondary is still a dwarf, despite it having nearly the 
same mass as the primary. 

Therefore, it is good to compare this lower limit for 
the radius with the radius that can be obtained from the 
photometric properties. According to Strassmeier et al. 
(pOOOT ), V=8 m ll, B-V=l™03, T Gff = 4790 K; and accord- 
ing to the Hipparcos catalogue (ESA 1997), ir = 4.74mas, 



0.9 km s _1 (Astronomical Almanach for 2002 (2001)) and 



corresponding to a distance of 211 pc. The V-magnitude 
could be too faint, if the star has been covered by spots at 
the time of the measurements and B-V could be too red; 
additionally, the interstellar extinction (Ay, E-q—v) is un- 
known. This means, that the resulting absolute magnitude 
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M v = l m 49 is an upper limit. With BC=CF43 (read off 
from Fig. 10.16, p. 197, in Gray 1992), the upper limit for 
the bolometric absolute magnitude becomes M DO i = 1™03, 
and the lower limit for the luminosity 29 Lq, leading to a 
lower limit of 7.85 Rq for the radius of the primary. This 
value is consistent with the lower limit from Eq. (^|) . Based 
on the good subtraction of the (3 Gem spectrum from that 
of HD 553, we could repeat the calculation using the tem- 
parature of (3 Gem: T cff = 4865 K (Drake & Smith |1991| ); 
this leads to (7.61 Rq), identical to the lower limit in Eq. 

While both results are lower limits, it would be quite 
a coincidence, if the difference of the rotation and orbital 
periods and the projection factor sin i on the one hand and 
the spot coverage and interstellar absorption on the other 
lead to the same change in the radius. It is more likely, 
that sin i Iot ss sin i orb w 1, P ro t ~ -Port,, V«V unspo ttcd, 
and Ay ss (despite a galactic latitude of just 2°). 

The effective radius of the Roche-lobe of the primary, 
i?RL, i.e. the radius of the sphere having the same volume 
as the Roche-lobe, is another quantity with which to com- 
pare the lower limit for the primary's radius. According 
to Eggleton QL983Q 



i?RL sin i — a sin i 



0.49 Q 2/3 



0.6Q 2 /3 + m (i + Qi/3) 
(8.859 ± 0.034) Rq 



(3) 



with the inverse mass ratio Q — mi/ma = 1.0429±0.0072. 

The fact that the Roche-lobe radius is only about 1 Rq 
larger than the radius of the primary shows that the as- 
sumption of synchronisation cannot be grossly wrong: If 
the rotation period were longer than lO^ (in the restframe 
of the primary) , the star would overfill its Roche-lobe and 
most likely form a shell around the whole system, which 
seems not to be the case. Thus, the deviation from syn- 
chronicity cannot be larger than 16%. 

If we consider synchronous rotation as valid, the star 
fills a considerable fraction (63.3% of the volume) of its 
Roche-lobe. Still, since the star is nearly 40 % smaller than 
the Roche lobe, we would not expect at present a large 
deviation from sphericity (see also Sect. [| where it was 
found that v sin i does not depend on phase). According 
to Gray (1992, App. B), a K0III star has on average a ra- 
dius of fl(KOIII) =11R®. Either the primary in HD 553 
is still expanding, in which case it will soon (in evolution- 
ary terms) fill its Roche-lobe and start losing mass towards 
the secondary, or its inclination makes R sin i — 7.6 Rq 
so much smaller than the typical radius, in which case 
the inclination would be about 44°. Since the system is 
eclipsing, the latter possibility seems highly unlikely; given 
the similar masses of the two components, but the very 
different luminosities, the former possibility seems to be 
very plausible. Thus, we expect HD 553 to become soon a 
semidetached system. 

Finally, the small errors in the radial velocity curve 
presented here, based on the long time span of 2 years, 
should allow a good transformation of future spectra into 



the restframe of the primary for surface imaging. However, 
many RS CVn stars show persistent active longitudes 
with activity cycles lasting several years (Berdyugina & 
Tuominen 1998). Thus, the relatively small amount of 
statistics, based on only 20 measurements over the com- 
paratively short time span of 2 years may not be enough 
to ensure that all systematic effects from spots on the ra- 
dial velocity curve are removed yet. Future spectra will 
therefore also be used to further improve the orbital pa- 
rameters. 
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